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'— '■ Abstract. The subject of this paper is to derive the solution of general- 

— , . ized fractional kinetic equations. The results are obtained in a compact form 

containing the Mittag-Lefner function, which naturally occurs whenever one 
is dealing with fractional integral equations. The results derived in this paper 
provide an extension of a result given by Haubold and Mathai in a recent 
q ■ paper (Haubold and Mathai, 2000). 

(N : 

O; 

^ ■ 1 Introduction and Preliminaries 

^— > ■ 

In terms of Pochammer's symbol 

> : 
X: 

^ ■ we can express the binomial series as 



( \ _ ( l,n=0 

\Ct)n — \ a (a+l)...{a+n-l),neN 



r=0 

The Mittag-Leffler function is defined by 



z n 
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This function was denned and studied by Mittag-Lefner (Mittag-Lefner, 1902, 
1905). We note that this function is a direct generalization of an exponential 
function, since 

Ei(z) := exp(z). 

It also includes the error functions and other related functions, for we have 
E l/2 {±z 1 ' 2 ) = e z [l + erf(±z 1 / 2 )] = e 2 er/c(^ 1/2 ), (3) 

where 

2 r z _ 2 

erf(z) ■■=- m \ e u du, erfc(z) := 1 - erf(z), z G C. (4) 

TT L ' Z JO 

The equation 

oo n 

gives a generalization of the Mittag-Lefner function. This generalization was 
studied by Wiman (1905), Agarwal (1953), Humbert (1953) and Humbert 
and Agarwal (1953) and several others. When (3 = 1, (5) reduces to (2). Both 
the functions defined by (2) and (5) are entire functions of order 1/a and 
type 1. A detailed account of these functions is available from the monograph 
of Erdelyi et al. (1955). The Laplace transform of E a ^(z) follows from the 
integral 

/ e- pt tP- 1 E a , (3 (\at a )dt = p- (3 (l-ap- a )- 1 , (6) 

J 

where Re(p) > |a| 1//Q , Re(f3) > 0, which can be established by means of the 
Laplace integral 

roc 

/ e-^-'dt = r(p)/pr, (7) 

Jo 

where Re{p) > 0,Re(p) > 0. The Riemann-Liouville operator of fractional 
integration is defined as 

a Drf(t) = ^ ^ f(u)(t - Uf-'du, V > 0, (8) 

with a D° t f(t) = f(t) (Oldham and Spanier, 1974; Miller and Ross, 1993; 
Srivastava and Saxena, 2001). By integrating the standard kinetic equation 

^N i (t) = -c i N i (t),(c i >0), (9) 
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it is derived that (Haubold and Mathai, 2000) 

N i (t)-N = -c i oD-'N^t), (10) 

where qD^ 1 is the standard Riemann integral operator. Here we recall that, 
in the original paper of Haubold and Mathai (2000), the number density of 
species i, iVj = iVj(t), is a function of time and Ni(t = 0) = N is the number 
density of species % at time t = 0. By dropping the index % in (10), the 
solution of its generalized form 

N(t)-N = -<? Dr"N(t), (11) 

is obtained (Haubold and Mathai, 2000) as 

^^\StSt?- (12) 

By virtue of (2) we can rewrite (12) in terms of the Mittag-LefHer function 
in a compact form as 

N(t) = N E u (-c u t u ), u>0. (13) 

In this paper we investigate the solutions of three generalized forms of (11). 
The results are obtained in a compact form in terms of the generalized 
Mittag-Leffler function and given in the form of three theorems. 

2 Generalized Fractional Kinetic Equations 

Theorem 1. If v > 0,/x > 0, then the solution of the integral equation 

N(t) - Not^ 1 = -c v Dr v N(t), (14) 

is given by 

N(t) = N T(fi)t»- l E^(-c u n, (15) 

where E Uitl (.) is the generalized Mittag-Leffler function defined by (5). 
Proof. We know that (Erdelyi et al., 1954) the Laplace transform of the 
Riemann-Liouville fractional integral is given by 

L{ Dt°f(t);p}=p-' T F(p), (16) 
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where 



F(p) 



-pu 



u=0 



f(u)du. 



Projecting the equation (14) to Laplace transform, it gives 

iv r(/i) 



N(p)=L{N(t);p} 



pf{l + (p/c)""}' 



By virtue of the relation 



L- 1 {p-p} = ^ y Re(p)>0, 



(17) 



(18) 



(19) 



it is found that 



N T(fi) 



p»{\ + (p/c)-"}_ 



NoT^L-Up^Y: 

{ r=0 



m-{p/c)- u r 

(r)! 



NoTirijti-iyc^L- 1 {p-^ r } 



r=0 

oo 



= N tt r^)t"-'E^-c"f). (20) 

The result (15) now readily follows by taking inverse Laplace transform of 
(18). For pi — 1, we obtain the result given by Haubold and Mathai (2000). 
Theorem 2. If v > 0, c> 0, d > 0, n > 0,Re(p) > c ^ d then for 

the solution of the equation 



N(t) - Not^E^i-trr) = -c v Dr v N(t), 
there holds the formula 



N(t) = ^o^—- [E„,»-„{-d v n - E^_ v {-c»n\ 



-v-\ 



(21) 



(22) 



Proof. Projecting (21) to Laplace transform and using (5) and (16), we obtain 

(-l)"(l) r (d/p)^ _ g (-l) r (l) r (c/p)" 



N p»-^ 
c v - d v 



.r=0 



(r)! 



r=0 



(r)! 
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Hence 

L-i{N(p)} = N(t) 



o 



c v - d v 
No 



Y j {-iyd rv L- 1 { P -^ +ri/ -^) - j^i-iyc^L- 1 ^- 

.r=0 r=0 

cc ( _ -^y ^jrvjji+rv—v—X oo ^ ^.ru-j-fi+ru— v— 1 



E 



E 



[E^-d^) - E v ,^ v {-c v t v )\ 



This completes the proof of (22). 
When fj, — u + l, theorem 2 reduces to 

Corollary 2.1. If v > 0, c> 0, d > 0, c ^ d, Re(p) > \d\ v / a , then for the 
solution of 

N(t) - N t v E v , v+1 {-d v t») = -c v D^N(t), (23) 
the following result holds 



N, 



N{t) = ^ Tu [E v {-d^)-E v {-c^)]. 



(24) 



On the other hand if d — > in (22), we arrive at 

Corollary 2.2. If c > 0, v > 0, fx > 0,Re(p) > \d\ v/a , then for the solution 
of 

V 



N(t) 



the following result holds 

N(t) 



= -c» DrN(t), 



E u ^ u {-c v n 



Note: When ji — v + 1, then for the solution of 



N(t) - N 



t v 



r(i/ + i) 



= -<? DrN(t), 



there holds the formula 



N(t) = ^ [1 - ^(-^Ol , 



(25) 



(26) 



(27) 



(28) 



where c > 0, v > 0. 

The case c = d is given by 

Theorem 3. If c > 0, v > 0, /x > 0, then for the solution of the equation 



N(t) - Not^E^i-JT) = -<? Dr, 
the following result holds 



(29) 



N (t) = ^V" 1 [E^-M") + (1 + v - ^E^i-c^)] . (30) 



Proof. Proceeding in a similar manner, it is observed that 

(-iy\2) r c ru L~ l { P -^- ru } 

r=0 



N(t) = N L- 1 {p-^l + c^r 2 } = NoYl 



(r)\ 



Not ^ r g ( _ 1)r [H(^ + ^-i) + (i + ,-^) }W 



iVr 



V 



oo / _„V4.v\r oo t _ n v+v\i 

V — ^ j , +(i + z/-/iV 1 c j 
^r(n/ + /i-i) 1 ^ 



= R.H.S. of (30). 



3 Conclusions 

The fractional kinetic equation (11) has been extended to generalized frac- 
tional equations (14), (21), (23), and (29). Their respective solutions are 
given in terms of the ordinary Mittag-Lefher function and their generaliza- 
tion, which can also be represented as FOX's H-functions. The ordinary and 
generalized Mittag-Leffler functions interpolate between a purely exponen- 
tial law and power-like behavior of phenomena governed by ordinary kinetic 
equations and their fractional counterparts, respectively (Lang, 1999; Hil- 
fer, 2000). A specific example for such behavior is the application of Tsallis 
statistics (Tsallis, 2002) to phenomena that may arise from fluctuations of 
temperature or energy dissipation rate (Lavagno and Quarati, 2002). The 
application of fractional kinetic equations to describe such phenomena has 
not been fully developed yet. 
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